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ABSTRACT 

We apply a new statistics, the factorial moment correlators, to density maps 
obtained from the APM survey. The resulting correlators are all proportional to the 
two point correlation function, substantially amplified, with an amplification nearly 
exponential with the total rank of the correlators. This confirms the validity of the 
hierarchical clustering assumption on the dynamic range examined, corresponding to 
0.5 h~"^Mpc — 50 h~"^Mpc in three dimensional space. The Kirkwood superposition with 
loop terms is strongly rejected. The structure coefficients of the hierarchy are also fitted. 
The high quality of the APM catalog enabled us to disentangle the various contributions 
from the power spectrum, small scale nonlinear clustering, and combinatorial effects, 
all of which affect the amplification of the correlators. These effects should appear in 
correlations of clusters in a similar fashion. 
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1. INTRODUCTION 

Correlation functions proved to be very useful statistics for studying the dis- 
tribution of galaxies because of their relatively simple connection to both theory and 
observations. The two-point correlation function gives only a rather limited view of 
the statistics, since it is equvalent to the power spectrum with all phase information 
lost. If the underlying distribution is non-Gaussian, the phases are correlated, and we 
need higher order correlation functions for the full description. In the case of galaxy 
catalogs the three and four point functions are quite strong (White 1979) displaying 
a clear sign of non-Gaussianity. Direct measurements of the higher order correlation 
functions have several difficulties: these quantities depend on an increasing number of 
variables (of which only a few are eliminated by symmetries), and are burdened with 
Poisson noise, biases, subtle edge and finite size effects, etc. Most of the high quality 
catalogs suitable for the determination of the higher order properties of the distribution 
are two-dimensional projections, and deprojection is only feasible under certain sim- 
plifying assumptions, while the redshift catalogs suffer distortions caused by peculiar 
motions. These are among the reasons why direct determination was only successful 
for the three-point funtion (Peebles & Groth 1975; Groth & Peebles 1977), and for the 
four-point function (Fry & Peebles 1978, Sharp, Bonometto & Lucchin 1984). Efforts 
were therefore refocused to indirect methods, of which the ones using moment corre- 
lators and/or moments were especially successful. Szapudi, Szalay & Boschan (1992, 
hereafter SSB) , and Meiksin, Szapudi & Szalay (1992) used both moments and mo- 
ment correlators to determine the higher order properties of the angular Lick and IRAS 
catalogs. Gaztanaga (1992), and Bouchet et al. (1993) used moments in cells to analyse 
the CfA and SSRS redshift surveys, and the IRAS galaxies with redshifts, respectively. 
All these references confirm the validity of a scaling law, the hierarchical assumption, 
defined in § 2.2. The three- and four-point functions from the Perseus-Pisces redshift 



survey, however, seem to support the inclusion of loop terms for a better fit (Bonometto 
et al. 1993). The moment methods generally determine an integral of the iV-point cor- 
relation functions over a cell (or a pair of cells, in the case of the moment correlators), 
and the results can be conveniently expressed in terms of the structure constants of 
the hierarchical ansatz, the Qjv's. The hierarchical ansatz was further supported by 
theoretical studies of the BBKGY equations (Davis and Peebles 1977, Fry 1984, 1986, 
Hamilton 1988), while it was also proposed on phenomenological grounds (Balian & 
Schaeffer 1989). The empirical constraints on the possible break-down of the hierarchi- 
cal ansatz are, however, quite uncertain, since redshift distortions in three dimensional 
surveys seem to show a similar trend (Lahav eA al. 1993, Vogeley eA al. 1994, Matsubara 
& Suto 1994), although other studies suggest that the effect is negligible within the 
measurment errors (Fry and Gaztanaga 1994). Projection could have a similar diluting 
effect, but according to our simulations, it is not the case, the Kirkwood assumption 
with loop terms could be easily distinguished from the tree hierarchy. The subject of 
the present study, the APM catalog, was also analysed by Gaztanaga (1994) using mo- 
ments in cells. He finds an excellent agreement with the hierarchy, and fits structure 
constants up to gg in a rather narrow range of angular scales 0.05° < 6 < 0.3°. For 
larger scales the uncertainties limit the number of points available to fit. We present 
a new method, based on the work in SSB, and the theoretical foundation in Szapudi 
& Szalay (1993a, hereafter SS), which can handle larger angular scales as well as small 
scales because of the nature of the factorial moment correlators, defined in the next 
section. This method has the additional advantage of automatically eliminating direct 
contributions from Poisson noise. Our aim is to fit for the hierarchy on angular scales 
0.1° — 5°, up to tenth order. The regular moments are incorporated in the normalized 
factorial moments, Wno, while new information is contained in the correlators Wnrm both 
defined in the next section. In §2. we define the statistics used, and briefly describe the 
method, in §3. we review the raw results obtained from the APM survey, §4. describes 



the fit for the hierarchy, in §5. we present the results of the fit, and the deprojection 
using all the maps, and in the final section we discuss and summarize the results. The 
Appendix contains some theoretical calculations which would have interrupted the flow 
of the main text, nevertheless they are needed for understanding and interpreting the 
results. 

2. DESCRIPTION OF THE METHOD 



2.1. Factorial Moment Correlators 
We introduce a set of novel statistics, the factorial moment correlators, deflned 



as 
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where we used the notatation {N)k = N(N — l)..(iV — 1 + A;) for the factorial moments 
of the counts in cells separated by a distance ri2- All these functions will go to zero 
when ri2 becomes much greater then the correlation length, since we subtracted the ap- 
propriate asymptotes. For a projected catalogue (as the APM survey, used throughout 
this paper) the distance between two cells is of course measured in angles. Similarly, 



we deflne the A;th factorial moments of the cell counts, denoted as w 
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where there is no subtraction. The statistic is a convenient reflnement of the (regular) 
moment correlators used by SSB with the addition that the factorial moments automati- 
cally eliminate the terms due to Poisson noise from the discreteness of the galaxy counts, 
since continuum moments become factorial moments after the transition to the discrete 
picture(SS). This use of factorial moments not only signiflcantly simplifles all the math- 
ematical formulae, rendering practical the study of much higher moments numerically 



than previously was possible, but also gives the statistic a very intuitive meaning, 

Wki ^ {pi...pip2---p2) - {p\){p[), (2.3) 

k I 

where p is underlying continuum random field (galaxy density) normalized such that 

its average is 1. The factorial moment correlator w^i is a special k -\- 1 order correlation 

function; the first cell is degenerate k times, the second is / times. Its two-point nature 

overcomes the difficulties of dealing with multidimensional quantities but the strength 

of the clustering signal is still preserved. Note, that the normalization follows from the 

continuum form. 

2.2. Analytical Model and Projection 

A useful property of the factorial moment correlators is that they can be expressed 
as a function of cell averages of the correlation functions up to a finite rank k-\-l. If P{x) 
is the generating function of the probability distribution, and P(x^y) is the bivariate 
generating function of the two-point distribution, the normalized factorial moments wok 
are generated (exponentially) by 

W(x) = P(- + l), (2.4) 

n 

where n = (iV), the average density, and the factorial moment correlators, Wki^ are 
generated (exponentially) by 

W{x, y) = P(^ + 1, ^ + 1) - P(^ + 1)P(^ + 1), (2.5) 

n n n n 

see SS for details. This general form could be used to fit for the integrals of the cor- 
relation functions over cells, but a further simplifying assumption can be taken. The 
hierarchical ansatz for the higher order correlation functions is defined as 

e(^) = Q^ J]e(ri)...e(riv-i). (2.6) 

trees 
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This is a Meyer expansion with tree level terms only. The symbolic summation means 
the following: we can imagine a tree graph in space using the N arguments of i^*^ -* as 
vertices. For each edge connecting the points ri^rj we multiply the N — 1 two-point 
correlation functions C (rij)- Finally we sum over all possible trees. Here only the 
functional form is fixed, and the Qn are structure constants. Note, that this assumption 
can be generalized to have different Qjv's associated with different graph types, or one 
can assume only the scaling described by the equation (Balian and Schaeffer 1989). By 
using the form of eq. (2.6) we essentially calculate combinatorially weighted averages 
for each N . In the next section we will show, that the hierarchical assumption is valid 
for the APM data in the dynamical range of our study. 

Under this ansatz the generating function takes the following form (SS): 

oo 

P{x) = exp Y, r^ix - ifQN, (2.7) 

JV=1 
where we introduced the following shorthand notation (SS) 

Tiv = ^^, ^ ^ , (2.8) 

where i^^ is the average of the two-point correlation in a cell (the same as i^ for the one 
point distribution), and v is the volume considered. If we denote the ('long') two-point 
correlation function between the two cells as i^; 

P{x,y) = P{x)P{y)[l + R{x,y)], (2.9) 

oo 

R{x, y) ~ 6 Y. (^ - l)'^(y - ^fQN+MTMrNNM. 

M=1,N=1 

Here R(x^ y) is approximate: only the terms linear in i^; were kept, since we assume that 
ii<:l. With this 

W{x, y) = P{- + 1)P(^ + l)R{- + 1, ^ + 1). (2.10) 

n n n n 



Expanding the generating function, the resuhing formulae agree with SSB if we set 
their Poisson-noise term to zero. This is a useful check since that result was obtained 
through an entirely different method. Since we operate with finite cells, the Qjv's in 
equations (2.7) and (2.9) are not exactly the same as in the hierarchical assumption. 
They involve correction factors, because we approximated the integrals of correlation 
functions (Balian and Schaeffer 1989, SS). The details of such a correction were worked 
out in Boschan, Szapudi and Szalay (1994). However, the correction factors in two 
dimensions are so small, that the measurement errors dominate, therefore it is not 
necessary include them in the calculations and in the fitting. 

Under the hierarchical assumption the higher order properties of the distribution 
are described by a set of Qjv's. Since the APM catalog is a two dimensional projection, 
we need to deproject the angular coefficients gjv to find the spatial Qjv's. This is possible 
via Limber's equation and the result is (Peebles 1980, SS) 

In = Jdlny{y'^{y)fy-''(''-'\ 

qN = QNl^-'lN/I^-'. (2.11) 

The APM luminosity function is the usual Schechter (1976) function, with variable M* 
and a, following Maddox et al. (1990a) 

M* = M*+Miz ; a = ao+aiz (2.12) 

M* = -19.8 M* = 1 ; ao = -1, ai = -2. 

The projection coefficients are found in Table 1. They are in good agreement with the 
calculations of Gaztanaga (1994), where he used a magnitude cut identical to one of 
ours. 



3. DATA — RAW RESULTS 

The APM Galaxy Survey has been described in detail by Maddox et al. (1990b, c). 
The survey is constructed from machine scans of 185 survey plates from the UK Schmidt 
telescope in a contiguous region which is relatively free from extinction, centered on the 
SGP. Plate regions containing bright stars, nearby galaxies, globular clusters or plate 
defects are masked off from the final data. Galaxy identification is automatically per- 
formed, giving a sample which is 90-95% complete with residual stellar contamination 
at the 5% level for 6j < 20.0, rising to 10% for 6j = 20.5. There is also residual con- 
tamination of the galaxy sample by merged images, which is estimated to be sa 5% in 
the magnitude range 17.0 < 6j < 20.0. 

We generated maps of the galaxy distribution in equal-area projections using 
pixel scales of 0.47°, 0.23°, 0.12°. Maps of the galaxy distribution in half-magnitude 
slices from 6j = 17.0 to 6j = 20.5, as well as two further maps covering 17.0 < 6j < 20.0 
and all galaxies brighter than 6j = 20.5 were considered for each pixel scale. A map 
of masked regions and survey edges was generated for each pixel scale to remove from 
our analysis all partially filled cells. Due to growth of empty regions as the resolution 
is decreased, the total area, and total number of galaxies available for analysis varies 
with pixel scale as given in Table 2., summarizing the most important properties of the 
different maps. For each density map we calculated the factorial moment correlators. 
Fig. 1. displays the results for the 0.47° cell size. The other cell sizes have similar 
appearance. Each figure displays the eight different magnitude cuts. The graphs show 
striking regularities. The parallel curves are the factorial moment correlators, Wki^ each 
corresponding to a degenerate k -\- /-point correlation function. 

There are three important features to notice on the graphs: (i) the parallelity of 
the lines, (m) the approximate equality of the moments corresponding to the same k-\-l^ 



regardless of the individual values of k and /, and, (Hi) the steady rise of the curves: i.e. 
there seem to be a linear amplification in log-space. This latter feature is even more 
apparent on Fig. 2. where only the (averaged) amplitudes are shown. 

3.1. Parallelity of the correlators 

The hierarchical assumption predicts parallel curves: there must be a long cor- 
relation function dominating the behaviour (see Eq. 2.9) Conversely, if the hierarchical 
assumption did not hold, loop terms would appear in the Meyer expansion of the higher 
order correlation functions, i.e. the approximation used to calculate R(x^y) would 
no longer be valid. The correlators could have terms with i^^^ (or even higher order) 
dominating at small scales, causing an upturn. The absence of this indicates that the 
hierarchical ansatz is a good approximation in the dynamic range considered. The pro- 
jection effects are subtle, however. Although a spatial hierarchy projects to an angular 
hierarchy, it is not true that any general functional form is preserved in projection. For 
instance let us consider the general Kirkwood approximation, 

if) 
F(^)(xi,...,XA.) = n^n(l + e^,), (3.1) 

where F^ ' is the disconnected correlation function. This contains loop terms, which 
appear as powers of ^ij higher than N — 1. In Appendix B we calculate the factorial 
moment correlator under the Kirkwood superposition. It exhibits the above mentioned 
behaviour: for a while the curves are parallel indicating the regime for which the Kirk- 
wood model is indistinguishable from the hierarchical assumption, and at small scales 
there is an upturn. This formula could be used to estimate the expected crossover 
scale in a spatial catalog, however, here we are dealing with an angular distribution. 
Although it is not possible to calculate the projection generally in this case, Toth et 
al. (1989) found that the cubic term, characteristic of the Kirkwood assumption for the 
three point function, will project to a similar term divided by the sum of the angles. 
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Since for large scales the Kirkwood superposition converges to hierarchical clustering 
with Qn = 1, there is a break-point where the cubic term starts to dominate the small 
scale behaviour, therefore we should see the upturn if our dynamic range includes the 
crossover scale. 

In the case of the factorial moment correlators the smallest scales in consideration 
arise when calculating the average over a cell. It could happen that while the hierarchical 
approximation is good for wi, it breaks down below the cell size. Then the parallelity 
of the curves would not be destroyed, because loop terms with w'j would be suppressed. 
Other loop terms, however, could affect the integral over one cell, introducing an effective 
Qjv, changing with the cell size. This effect could manifest itself in Q^ and higher 
terms, causing a systematic shift. There is no sign of this systematics in the APM data 
within the errors, meaning that we did not detect any loop turn down to the scales 
corresponding to the cell size. 

Since an analytic formula is available only for the cubic term, it is very hard to 
estimate the crossover scale, although we conjecture that it would be well within our 
dynamic range for the higher moments, where the deviation is expected to be the most 
significant. The equation for a possible numerical projection can be obtained from the 
generating function of w^i as expressed in terms of the disconnected correlation functions 
(see Appendix B): the factorial moment correlators are essentially the integrals of the 
disconnected correlation functions. After the substitution of the Kirkwood superpositon 
and, changing the measure dx — )• dQr'^drtp(r) = dji, where tp is the selection function, 
O is the solid angle of the pixel (see SS for details), we obtain the following equation 
(to be compared with Limber's equation) 



I dij.i . . . dij.k n (i+czj) (3.2) 

pairs 

/ (i/ii . . . dfj^k / dfj^k+i ■ ■ ■ djik+i n (1 + iij)- 
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These equations are increasingly divergent with increasing order, since e.g. for wno 
the leading term in the integrand is ~ r )'"("-i)/2 while the measure is ~ d^ r. The 
correlation function, however, cannot be a true power law, there has to be a small scale 
cut-off. This renders the integral convergent, but its value depends sensitively on this 
largely unknown cut-off. If we overestimate this cut-off, we can still get a lower bound 
on the integral. A direct Monte-Carlo approach with the Kirkwood approximation 
is highly unstable, since the most important contributions to the integral arise from 
very unlikely configurations. Therefore we developed a projection method based on 
Riemannian sums, described in Appendix D. We projected i^(r), the two point function, 
and 1^50 using both the hierarchy and the Kirkwood assumption with a certain cut-off. 
We have chosen 1^50, because among the factorial moments this has the highest degree 
of divergence in the integral, and the smallest available cell size (0.12°), where the 
most prominent effect is expected. We used the 17 — 20 magnitude slice of the catalog 
for these numerical studies. Fig. 3 shows the projection of the two-point function 
with and without a cut-off. Without smoothing, the resulting correlation function is 
highly consistent with the APM two-point function (Maddox eA al. 1990c). Clearly from 
Fig. 3., a cut-off of 0.5h~"^Mpc is rejected by the observed APM two-point correlation 
function. We therefore use this value for the cut-off upper limit to produce a lower 
bound on the integral. The Kirkwood assumption predicts Wko > 1820 (see Appendix 
D. for details). The observed value is ~ 14.7, therefore the APM data strongly reject the 
Kirkwood assumption. For a hierarchy with Qn = 1 we used the same code to project 
1^50 for a comparison. In this case there are no convergence problems, so we obtain an 
approximation, rather than a lower bound. The result is 7.6, while using the projection 
of Qn = 1 and the expression in Appendix C. gives 6.9. This is fully consistent with the 
hierarchy, since if we perturb the Qjv's slightly, we can easily obtain a value very close 
to the measurements. In summary, the factorial moment 1^50 fully supports a hierarchy 
with Qn — 1, and strongly rejects the Kirkwood assumption. 
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3.2 Approximate Degeneracy, Amplification 

If the degeneracy and the hnear amphfication in log space were exact we could 
conclude that the hierarchical structure constants satisfy gjv = 1 for all N (see Appendix 
A). These properties of the correlators are, however, only approximate. Even if Qn = 1 
were true for the spatial coefficients, projection would modify it through eq. (2.11), 
although not by orders of magnitude. Since we want to consider gjv = 1 as a 'ground 
state' let us define the following quantity: 

fun. = - ( "'"'" ) . (3.3) 

Wl \nmWnOWmO J 

Were the gjv's all equal to one, fnm = 1 exactly would hold through all order (see 
Szapudi & Szalay 1993b, hereafter SS2, or the expansion of Appendix C). Since we 
expect this to be only an approximation we compiled a plot of fnm for all of the density 
maps obtained from the APM survey. From Fig. 4. the deviation is less then a factor 
of two, so it is indeed a good approximation. According to this the correlators can be 
written in the following form 

Wnm = nmWlWnOWmO fnm- (3.4) 

Here the wi contains the information from the power spectrum, and that determines 
angular dependence of the correlators, nm is a purely combinatorial factor, WnoWmO 
depends on the gjv's up to max(?7z, n) and the Wg^ and contains the vast majority of the 
small scale nonlinear clustering. The corrective factor fnm is the residual accounting 
for fine details of the small scale clustering and is found to be of order unity. This is a 
remarkably good approximation if we consider that the amplification of the correlators 
ranges 3 — 4 orders of magnitude, and eq. (3.4) explains it up to factor of 2. This 
shows that the amplification process is very robust, it depends on the details of the 
clustering hierarchy only subtly. These considerations are valid for the amplification of 
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the cluster correlation function, since Wnn/w'^Q and w^i/wno correspond to the cluster- 
cluster and cluster-galaxy correlation functions in a catalog where cluster selection is 
modeled by a soft clipping, i.e the catalog is raised to a power (SSB). This is demon- 
strated on Fig. 5, where the different contributions to both the Wko^s and Wki^s are 
plotted: (iV)ytPjv/(iV) . The curves are increasingly shifted toward higher values of N 
indicating, that the most important contribution to the different moments comes from 
cells with higher and higher counts. Therefore raising the catalog to a power (or bet- 
ter replacing with a factorial moment) indeed corresponds to a contrast enhancement 
similar to biasing (Kaiser 1984, Bardeen et al. 1986). Eq. (3.4) predicts that the shape 
of the cluster correlation function is identical to the galaxy correlation function and 
the amplification is roughly n^ and n respectively, for the auto and cross correlation 
function, and that the amplification depends on the details of the selection algorithm, 
increasing with richness. These results are very robust since fnm is found to be of order 
unity, therefore calculations assuming gjv = 1 (SS2) yield very realistic results. An 
approximate degeneracy and linear amplification in log space can arise simply as a com- 
binatorial effect: if we plot k\l\ logarithmically we obtain similar graphs. (log(A;!/!) ~ 
A;(log A; - 1) + /(log / - 1) + log A;/2 + log //2 + 2 log ^Ivr ~ (A; + /)(1 + 0(log / + log k), 
via the Stirling formula). This means that Wko ~ (A; — 1)! according to eq. (3.4), which 
is a good approximation from Fig. 2. 

One should be aware, that the 'cluster correlation' as derived in this paper from 
the normalization of the factorial correlators is only an approximation to correlation 
functions based on more realistic cluster finding algorithms (e.g. see Dalton eA al. 1992). 
Here we use distinct, fixed cell sizes, use only projected counts, and discard detailed 
magnitude information. In our approach, a single cluster extending over several cells 
would be counted with multiple weight, thus a one-to-one comparison with an objective 
cluster finding algorithm is not trivial. On the other hand, our approach makes an 
analytic understanding of the phenomenon very easy, and the gross features i) identical 
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shape of i^cc-, Ccg, and i^^^, ii) i^cg — yCggCcc should be very robust and vahd for real 
cluster surveys as well. Our method shows, that the actual amplification has to depend 
on the details of the selection algorithm, like the cluster radius and threshold, and these 
will differ for every cluster catalog, but the above relation should persist. 

To summarize considerations in this section, in the examined regime 0.1° — 5° 
the data are well described by a hierarchy with the g^'s close to 1. Since only the 
hierarchy is preserved by projection, we performed simulations, if projections could 
result in a similar behaviour in an angular catalog, even if the spatial disribution is 
not hierarchical. The most motivated and well defined hypothesis with loop terms, the 
Kirkwood superposition can be rejected by these data. The observed amplification of 
the correlators is disentangled into factors due to the power spectrum, the small scale 
non-linear clustering, combinatorial effects. 

4. FIT 

To quantify the higher order properties of the APM catalog we used the raw 
factorial moment correlators to fit for the gjv projected hierarchical structure constants. 
We performed a maximum likelihood fit similarly to SSB using the hierarchical model 
to calculate the Wki^s for a given set of gjv's. We assumed a Cauchy distribution for 
the errors because of its robust handling of the outlier points (Press et al 1987). We 
calculated all the Wki^s and wqi^s for a given set of g„, and a set of i^^'s for the catalogs 
which had the same magnitude cut but different cell sizes. This was compared to the 
measured factorial moment correlators from the density maps, the likelihood function 
was calulated, and the g^'s and i^^'s were adjusted by the AMOEBA algorithm (Press 
et al 1987). The weights for each term in the likelihood function were assigned using 
a particular set of g^'s, and at end of the fit we checked if the resulting g^'s settled 
reasonably close to the ones assumed for the weights. We performed simultanious fits 
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for a given magnitude cut with all the different cell sizes. Then we ran AMOEBA 
several times with random initial conditions, to assess a formal error on each gjy. The 
results are presented in Table 3. The coefficients of Table 1. can be used to convert 
these to spatial coefficients. The lower order q^ are much better determined than the 
higher order ones, because the lower order ones occur in more terms (see appendix C.) 
then the higher order ones, therefore the algorithm is more sensitive to changing those 
parameters. 

5. SUMMARY AND DISCUSSIONS 

We used density maps obtained from the APM galaxy catalog to calculate facto- 
rial moment correlators. The resulting graphs showed remarkable regularity over a wide 
range of cell sizes and magnitude cuts: all the curves are parallel to the two point corre- 
lation function and there is a uniform amplification in log space logWki — const(A; -|- /). 
The former points toward the validity of the hierarchical assumption for the dynamic 
range considered, and latter is supported by both by the gjv's being not far from 1 and 
by combinatorial conspiracy. We fitted for the structure constants gjv, and used the 
APM luminosity function to deproject them for the different magnitude cuts. Table 3. 
contains the results for the different magnitude cuts, the errors cti are calculated for- 
mally from hundred fits from random initial conditions. The dispersion was calculated 
cumulatively after each fit. We found that the errorbars and the average are sufficiently 
settled, therefore further runs with random initial conditions would not have given more 
information or tighter errorbars. In an ideal case the different magnitude cuts would 
give the same Qjv's when deprojected: this gives a chance to give a better estimation 
of the errors. Ei is the deprojected average random error, and E2 is the dispersion 
calculated from the different magnitude cuts. This latter error might contain several 
contributions: systematic errors from the possible break down of the hierarchy, different 
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amount of contamination for the different density maps, possible errors in the luminos- 
ity function (especially at the faint end). In general any qualitity or error of the data 
that changes with the depth can influence the flnal flt. Although these sources of error 
cannot be separated now, Ei + E2 gives probably the most realistic error estimation up 
to now. 

We conclude that we found the validity of the hierarchical assumption over a 
scale range of 0.1° — 5°. For the same range we fltted the structure coefficients which 
are remarkably close to 1, although the higher coefficients are also consistent with zero. 
This result confirms the validity of the hierarchy on scales 0.5h~"^Mpc — 50h~"^Mpc, 
while it is inconsistant with the Kirkwood superposition on the same scales. 

APPENDIX A: AN ANALYTICAL MODEL 

Let us assume that the amplification is exactly linear and the degeneracy is per- 
fect, and pursue the consequences. These conditions are equivalent of Wki in the form 
^i{A/n) ~^ . /,From this the generating function of Wki can be evaluated as 

^^%^ = (e^^ - l){e^y - 1), (A.l) 

yielding a very strong constraint for the structure constants: they have to be power law, 

Qn = Qs~\ N>2. (A.2) 

This can be proved by induction using the separability of the generating function. This 
would mean essentially the very attractive assumption Qn = 1, since Qi = Q2 = ^■ 

If we relax the condition of linear amplification, but keep the degeneracy of A; -|- /, 
the generating function becomes of the following form 

W{x,y)=B{x+y)-B{x)-B{y) (A.3) 
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where B(x) is the exponential generating function of b(k + /) = Wki- Note, that in the 
former case B{x) = ^i{e — 1). Of course the conditions are not exact, therefore more 
general forms are allowed, but nevertheless this gives a very useful toy model (see SS2 
for calculation of the cluster correlations within the framework of this model), and a 
rough idea of the order of the Qjv's. 

APPENDIX B: FACTORIAL MOMENT CORRELATORS: KIRKWOOD 

To calculate the factorial moment correlators we need to express the generating 
functions with disconnected correlation functions, since the Kirkwood superposition 
is natural in those terms. The path described here briefly is analgous to the methods 
developed in SS, where connected correlation functions were considered. The generating 
functional expands in terms of the reducible correlation functions F^ ' (no exponential), 

Z[J] = J2 -^jlr / ^^1 • • • ^^^^^^"^^(^1 , • • • , xn)J{xi). . . J{xn). (B.l) 

Jv=o ■ "^ 

The generating function, similarly to SS can be expressed as P{z) = (e''-^^^~^'), where 
^r{x) = J e{y)WR{x,y)dy is the flltered continuum fleld, Wr is the fllter function, 
and e(y) is the underlying random fleld. If we take the source function J*(x') = 
Wr(x^x')(z — l)/z, and Wr is the characteristic function of the volume considered, 

^(^)= T. m dx,...dxNF^^\x,,...,XN). (B.2) 

JV=1 "^^ 

Similarly the bivariate generating function, P{zi,Z2) = (^e''-^^^^~^'''-^^^^~^') , can be ex- 
pressed as 

{Z, - 1)^{Z2 - 1)^ 



AT A/f ^ 



N\M\ 

N,M 

dxi . . . dxN / dxN^i . . . dxN^MF "*" (^^i , • • • , xn^m) { ■ (B.3) 
Vi Jv2 



17 



According to eqs. (2.4) and (2.5) the integrals in the previous two equations are n Wko = 
{{N)k) and n ~^ Wki = {{N)k{N)i). The previous two equations together with eqs. 
(2.4) and (2.5) mean that the factorial moments and the factorial moments correlators 
measure the avereges of the disconnected correlation functions over a volume or two 
disjoint volumes respectively. 

If we substitute the Kirkwood superposition of eq. (3.1) in the first integral 

. (^) 

{{N)k) = n^ dx^... dxk 1[{1 + U)- (B.4) 

i=l 

The second integral similarly yields 



{{NMN)i)=n 



k+i 



/ dxi...dxk / dxk+i ...dxk+i TT (1 + Cafe)- (B-5) 

JVr JV2 ,^1 



These equations are needed when calulating the projection of the factorial moments 
correlators. 

APPENDIX C: FACTORIAL MOMENT CORRELATORS: HIERARCHY 

Under the hierarchical assumption the Wki^s can be expressed finitely as a function 
of the average of the correlation function i^^, the 'long' correlation function i^; and the 
structure constants q^- The Wno^s are obtained from the expansion of the exponential 
generating function eq. (2.4) as the coefficient of x"" /n\: 



"^20 = 1 + 6 

W40 = 1 + 6 ,^5 + 3 ,^5^ + 12 ^s'^ q3 + 16 C/ ?4 

u;5o = 1 + 10 6 + 15 6' + 30 6' 53 +30 6' 53 + 80 6' 54 + 125 e/ 55. 



c 3 



The Wnk are obtained in a similar fashion through eq. (2.5): 

Wii/wi = "i + "i is + Q is qi + '^ is^ qA 

Wiijwi = 4 + 12 e. + 12 6 53 + 24 is" qz + 36 6' qi + 64 C. 
w^5i/w^; = 5 + 30 6 + 15 6' + 20 6 gs + 120 6' gs + 60 6' ^3' + 

90 6' 54 + 170 e/ 54 + 320 6' 55 + 625 C/ ge 
W22IW1 = 4 + 8 ^5 53 + 4 ^^^ 54 
W32/W1 = 6 + 6is + 'i-8isq3 + 

6 6^ ?3 + 30 ^^^ 54 + 18 C/ qs 
Wi2/wi = 8 + 24e. + 32 e. 53 + 72 C.' 53 + 24 C/ 53' + 

96 6' 54 + 24 e/ 54 + 200 e/ 55 + 128 C/ 56 
u;52/w^; = 10 + 60 6 + 30 6' + 50 6 53 + 300 6' 53 + 30 6' 53 + 

240 e/ 53' + 220 6' 54 + 460 6' 54 + 280 C/ 53 ^4 + 

820 e/ 55 + 180 e/ 55 + 1890 e/ 56 + 1250 C/ 57 
«^33/w^; = 9 + 18 6 + 9 6' + 36 6 53 + 36 C/ 53 + 90 C/ 54 + 

54 e/ 54 + 108 e/ 55 + 81 e/ 56 

u;43/w^; = 12 + 48 6 + 36 C/ + 60 6 53 + 180 C/ 53 + 
72 e/ 53 + 72 e/ 53' + 216 6' 54 + 288 6' 54 + 

108 e/ 53 54 + 516 e/ 55 + 300 e/ 55 + 708 e/ 56 + 576 e/ 57 

u;53/w^; = 15 + 105 6 + 135 6' + 45 6' + 90 6 53 + 600 6' 53 + 
450 e/ 53 + 540 e/ 53' + 180 C/ 53' + 435 6' 54 + 
1410 e/ 54 + 645 e/ 54 + 1380 C/ 53 ^4 + 720 C/ 54' + 
1680 e/ 55 + 2040 e/ 55 + 540 C/ 53 ^5 + 4605 C/ 56 + 
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2685 e/ ge + 6630 C/ qi + 5625 C/ gs 

u;44/w^; = 16 + 96 6 + 144 C.' + 96 6 gs + 480 6' gs + 576 C/ ^3 + 
288 e/ gs' + 432 e/ ga' + 432 6' ^4 + 1152 C/ ^4 + 
144 e/ 54 + 864 e/ 53 54 + 1376 6' ^5 + 2400 C/ gs + 
1536 e/ 53 55 + 2832 C/ ge + 1536 C/ ge + 4608 C/ ^7 + 4096 C/ gs 

u;54/w^; = 20 + 180 6 + 420 6' + 180 tJ" + 140 6 ^3 + 1200 6' ^3 + 
2340 e/ 53 + 360 e/ 53 + 1200 6' ^3' + 2880 C/ ^3' + 
720 e/ 53' + 780 6' 54 + 3800 6' qi + 3300 C/ ^4 + 
4920 e/ 53 54 + 3720 C/ ^3 ^4 + 2880 C/ ^4' + 3400 6' ^5 + 
10080 e/ 55 + 2040 e/ 55 + 9840 C/ ^3 ^5 + 5120 C/ ^4 ^5 + 
10860 e/ ge + 18420 C/ ge + 11340 C/ ^3 q& + 24780 C/ ^7 + 
13260 e/ 57 + 42980 C/ gs + 40000 C/ go 

u;55/w^; = 25 + 300 6 + 1050 6' + 900 C/ + 225 C/ + 200 6 
2400 6' 53 + 7800 e/ 53 + 3600 C/ ^3 + 3000 6' ^3' + 
14400 e/ 53' + 1800 e/ 53' + 7200 C/ ^3' + 1300 6' ^4 + 
9500 e/ 54 + 16500 is" qi + 5100 C/ ^4 + 16400 C/ ^3 ^4 + 
37200 e/ 53 54 + 13200 C/ ^3' ^4 + 
14400 C/ 54^ + 20800 C/ 54^ + 6800 6 
33600 C/ 55 + 20400 C/ gs + 49200 C/ 
10800 e/ 53 55 + 51200 e/ 54 55 + 27150 C/ ^6 + 
92100 e/ ge + 26850 C/ ge + 113400 C/ ^3 ^6 + 
100000 e/ 54 ge + 82600 C/ ^7 + 132600 C/ ^7 + 75000 C/ ^3 ^7 + 
214900 C/ gs + 112500 C/ qs + 400000 C/ go + 390625 C/ gio 



*g3 + 

' ^ I 

53 + 



^ ^3 ^5 + 
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These formulae compare directly with the formulae listed in SSB without Poisson 
noise, 1/iV = 0. The simplification achieved by the factorial moments can be illustrated: 
the W42 term contains 23 terms if one uses regular moments as opposed to 9 for the 
factorial moments. The higher we go, the more dramatic is the difference due to the 
combinatorial explosion of terms. 

APPENDIX D: PROJECTION UNDER THE KiRKWOOD ASSUMPTION 

As described in the main text the integration of Eqs. (3.2) for the projection 
would be highly unstable with a Monte Carlo approach. We therefore used a direct 
Riemannian sum to approximate the integral. The integration volume is a pyramid 
corresponding to the solid angle O. We have chosen the smallest cell size (0.12°) and 
the 17 — 20 magnitude slice for practical investigations. We divided this 'pencil beam' 
into 2 slices of equal length, and on each slice we approximated the correlation function 
with its average on the slice. In principle, of course, we could have used any value of 
the integrand within the cell, but we anticipated that the average might give a more 
stable performance. Although the method is completely general, we used it to calculate 
W50 (where the effect is expected to be the most prominent), and we describe it in those 
terms for simplicity. The generalization should be straightforward, although tedious 
for higher order factorial moments, or moment correlators. In the latter case the only 
difference is, that the integration volume comprises of two pencil beams. 

The projection of Wko is a five fold integral. The five variables are distributed 

among the cells. There are 7 distinct cases, according to how many variables fall within 

one cell: 

5 

4 1 

3 2 

3 1 1 

2 2 1 

2 111 

21 



11111. 

These are all the possible partitions of 5 into integers. 4 1 means that four points are 
in one cell, and one point is in a different cell, etc. The sum of the contributions over 
all possible configurations with the appropriate measure calculated from the selection 
function will approximate the integral as the cell size approaches zero. The problem 
with the Monte Carlo approach lies in the fact that the most important contribution to 
the integral comes from the first partition 5, where all the points lie in the same cell. 

The average of the correlation function over a cell was calculated by a Monte Carlo 
simulation, and the following formula was found to be an adequate approximation, 

where the pyramid slice was approximated with a rectangular parallepiped with the 
ratio of the short and long side x = d/l, and we used 7 = 1.7 and vq = 5h~"^Mpc. 

The following notation is needed for the approximation of the correlation function 
1^(1 A^' — A^' I) between two different cells (the upper index in parentheses numbers the 
vector, while lower index defines the coordinate). Let us draw the vector x^^' and x^^' 
to the center of mass of the two cell. In a general calculation these could be in two 
different pencil beams. Then we can write 

^ = ^(1) _ ^(2) ^ ^(1) _ ^(2) ^ ^(1) _ ^(2) p_2) 

where s^^' points from the center of mass, so its average {s^^') = 0. The Taylor expansion 

a\ r I) = aR) +s^^ \s=o +ls.s,^^ |.=o + . . . , (D.3) 

OSi 2 OSiOSj 

where R =| x^^' — x^^' \ the distance between the center of mass of the two cells, 
X = x^^' — x^^K s = s^^' — s^^K and the Einstein convention is used. After differentiation 



a\r\) = aR) + -^ + ^s.s, 






1^(^"-i^')+^4 



(D.4) 
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If both of the ceUs are in the same beam, we can calculate the average of this using 
{s,sj) =6,j{s1) as 



{i) = im + \{s]) 



x^ 1 f 



(D.5) 



If we substitute R^' = —j^ and R^^" — R^' = 7(7 + 2)i^, and we use that the two cells 
are in the same pencil beam therefore x = (0, 0, i?), 

(0 =e(i?) - 1 {{{s^'^f) + {{s^'^f)) ie{R)/R+ (D.6) 

This correction falls off as l/i?, and is therefore important only for the few nearest 
neighbours. For the touching cells, however, it is not accurate enough, and we used the 
constraint for the total correlation of the two cells. If xn and 2:22 denote the averages 
of the correlation functions in the two touching cells, a;i2 is the average of the cross 
correlations between them, and 1^0 is the average of the total correlations in the double 
cell, we have the following formula because of the additivity of the integral: 

4^0=61+62 + 262. (D.7) 

In the present case the two cells are nearly identical, because of the small angle of the 
cone, therefore 61 — 62- 

Now to calculate the integral 1^50 = J djii . . .dji^ Ylpairsi'^ + 6j) '^6 needed to 
calculate the correlations on the set of cells, and the discrete integration measure pi 
replacing dfx was simply calculated on each cell using the selection function such that 
^Pi = 1 hold. When implementing it as a computer program we need only as many 
loops as the number of integers 5 was partitioned into, e.g. we need only one loop for 
the partition with all points in the same cell. If the integral is finite, the approximate 
sum is convergent as the cell size becomes smaller. Similarly the hierarchical W50 could 
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be projected using a formula similar to Appendix C, except we left gi, 52 in the formula 
explicitly, although they both equal 1, they were needed to locate the appropriate 
orders. Then it was a simple matter to assign the correct power of the measure and the 
correlation function for each term. For instance 30g2<Z3Cs '^^^ associated with ^Opjp^j^i^^ , 
where i and j denote the loop variables, etc. Finally, the same code calculated the 
projection of the two point function and calculated Wg for mainly test purposes. 

The results of the integration are the following. The two-point function and the 
average of the correlation function over the APM cell are well recovered by the code. 
The projection of W50 under the hierarchical assumption with all Qn = 1 converged at 
cell sizes as big as 96h~"^Mpc! The integration gave 7.6, in reasonable agreement with 
the formula in Appendix C. combined with the projection of the Qjv's from Table 2. 
gave 6.92. The measured value is 14.7. The factor of two difference can be accounted 
for by letting the Qjv's differ from 1 by a small amount. For the Kirkwood assumption 
we can only expect convergence if we use a cut-off in the correlation function, otherwise 
the integral is divergent. Every approximate sum until convergence is therefore a lower 
bound, and if we use a larger than real cut-off we obtain an even lower bound. The 
cut-off from the APM correlation function has to be smaller then 0.5h~"^Mpc. Using 
that as a cut-off and calculating the integral with smaller and smaller cell sizes, we 
did not reach convergence up to 24h~"^Mpc cell size, but the lower bound obtained is 
1820, some 124 times the observed value. Therefore the data safely reject the Kirkwood 
assumption. 
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Table Captions 

Table 1. — For each magnitude cut {nimin — i^max) and for three different 
ceU sizes (/) the average number of galaxies in a ceU (Nave) is displayed, ^^st is the 
estimated amplitude of the correlation function using only second order quantities, ^fn 
is the result of the fit up tenth order. The fact that these two determination of i^ do not 
differ significantly supports the validity of the fit. For the sake of completeness we give 
the total number of objects (Ntot)} and the surface area (A) covered by each catalog 
with different cell sizes. 

Table 2. — The coefficients used for deprojection of the angular structure con- 
stants to spatial ones. For illustration one of the magnitude cuts is calculated with 
both h = 1 (the dimensionless Hubble constant) and h = 0.5, the rest of the calculation 
h = 0.5 is used since the results are not sensitive to this parameter. 

Table 3. — The main numerical results of the paper are displayed. The fitted 
angular coefficients q^ are shown for each magnitude cut, the formal error from hundred 
runs with random initial conditions is cti. The spatial Qjv's are calculated using all the 
angular gjv's and the coefficients of Table 2. Ei is the formal deprojected error, E2 is 
the dispersion from the different deprojections used when calculating the average. 



Figure Captions 

Fig. 1. — The factorial moment correlators from the catalog of 0.47° cell size 
are shown for each magnitude cut. The magnitude cut is coded under the figures in a 
way that e.g. 170200 means a cut between 17.0 and 20.0. The graphs are strikingly 
parallel. The correlators from the other density maps with cell sizes 0.23°, and 0.12° 
have similar appearance. 

Fig. 2. — The amplitude log(wki/wii) is shown as a function of A; + /. There is a 
small splitting in the degenaracy in A; + / compared to the average trend which is nearly 
linear in log-space. The different cellsizes are coded as follows: square: 0.12°, triangle: 
0.23°, and cross: 0.47°. The solid line corresponds to the normalized factorial moments 

Fig. 3. — The projection of the two-point correlation function is shown with 
a smoothing of a = 0.5h~"^Mpc and without smoothing. The projected curve with- 
out smoothing is consistent with the APM two-point function, while the other curve 
is strongly inconsistent. Therefore the cut-off in the real correlation function is less 
then 0.5h~"^Mpc. The w^q moment, calculated from the Kirkwood assumption with 
0.5h~"^Mpc cut-off, is orders of magnitude larger then the observed value (which is 
consistent with the hierarchy). The Kirkwood assumption is strongly rejected, because 
smaller cut-off would result in even larger w^q. 

Fig. 4. — The quantity log(fnm) is displayed as a function of n -\- m. This plot 
enables us to disentangle the different effects contributing to the amplification of the 
correlations. Note that fnm is of order unity while the amplifications on Fig. 2. are of 
order 10^ ^10^ 

Fig. 5. — log((N)kPN/{N)'') versus N is plotted for A; = 0, . . . , 5 to show that 
the factorial moments correspond to a contrast enhancement. For this plot the 17.0 



to 20.0 magnitude cut with 0.23 cell size was used. The most important contribution 
indeed shifts towards higher iV's as k increases. 
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Table 1. 



Density Maps from the APM Survey 



/ 




0.117188° 






0.234375° 






0.468750° 




AT 




2264402 






2182470 






1981523 




Atot 




4227 






4069 






3692 




f^min f^max 


^ ' ave 


iest 


ifri 


A^ 


iest 


if^i 


A^ 


iest 


if^i 


00-20.5 


6.791 


0.165 


0.181 


27.192 


0.106 


0.122 


108.88 


0.0665 


0.0665 


17.0-20.0 


3.986 


0.224 


0.224 


15.961 


0.147 


0.158 


63.92 


0.0948 


0.0949 


17.5-18.0 


0.158 


0.705 


0.703 


0.633 


0.466 


0.519 


2.539 


0.295 


0.295 


18.0-18.5 


0.308 


0.504 


0.504 


1.235 


0.333 


0.424 


4.953 


0.224 


0.252 


18.5-19.0 


0.572 


0.338 


0.338 


2.290 


0.226 


0.211 


9.165 


0.147 


0.167 


19.0-19.5 


1.072 


0.234 


0.229 


4.292 


0.154 


0.194 


17.20 


0.101 


0.105 


19.5-20.0 


1.798 


0.166 


0.166 


7.198 


0.107 


0.115 


28.82 


0.0678 


0.0681 


20.0-20.5 


2.713 


0.122 


0.122 


10.863 


0.0755 


0.0757 


43.491 


0.0460 


0.0459 



Table 2. 



Correction Factors for Tree Graphs up to iV = 8 



i/7 


1.5 


2D 

1.65 


1.8 


1.5 


3D 

1.65 


1.8 


1 


1.000 


1.000 


1.000 


1.000 


1.000 


1.000 



Table 2. 



Projection Coefficients for the Qn's 



f^min 


f^max 


h 


i?3 


ii4 


i?5 


i?6 


i?7 


Rs 


-Rg 


-Rio 


17.0- 


-20.0 


1.0 


1.161 


1.434 


1.830, 


2.386 


3.157 


4.221 


5.690 


7.719 


17.0- 


-20.0 


0.5 


1.178 


1.483 


1.934 


2.579 


3.492 


4.781 


6.601 


9.174 


00.0- 


-20.5 


0.5 


1.168 


1.453 


1.869 


2.454 


3.271 


4.405 


5.981 


8.171 


17.5- 


-18.0 


0.5 


1.159 


1.426 


1.813 


2.351 


3.093 


4.110, 


5.505 


7.420 


18.0- 


-18.5 


0.5 


1.157 


1.419 


1.797 


2.323 


3.045 


4.033 


5.382 


7.228 


18.5- 


-19.0 


0.5 


1.153 


1.411 


1.780 


2.293 


2.994 


3.950 


5.251 


7.025 


19.0- 


-19.5 


0.5 


1.151 


1.402 


1.763 


2.261 


2.941 


3.863 


5.115 


6.813 


19.5- 


-20.0 


0.5 


1.147 


1.393 


1.744 


2.228 


2.885 


3.775 


4.976 


6.600 


20.0- 


-20.5 


0.5 


1.144 


1.384 


1.726 


2.196 


2.831 


3.688 


4.840 


6.391 


17.0- 


-17.5 


0.5 


1.162 


1.433 


1.827 


2.377 


3.136 


4.181 


5.618 


7.596 



Table 3. 



Results of the Fits for the g^'s 



ITT-min 


ITT-max 


order: 


3 


4 


5 


6 


7 


8 


9 


10 


17.0- 


-20.0 


qN 


1.160 


1.964 


5.274 


9.900 


2.123 


1.236 


-0.653 


0.380 






cri 


0.002 


0.090 


0.625 


2.459 


9.208 


9.053 


9.728 


11.452 


00.0- 


-20.5 


qN 


1.141 


2.224 


6.812 


2.963 


1.447 


0.971 


0.066 


0.868 






cri 


0.005 


0.118 


1.110 


4.363 


5.254 


5.524 


5.486 


5.479 


17.5- 


-18.0 


qN 


1.624 


3.056 


5.868 


12.192 


6.035 


3.118 


-0.197 


3.847 






c^i 


0.002 


0.009 


0.124 


1.215 


14.664 


17.136 


21.994 


20.197 


18.0- 


-18.5 


qN 


1.373 


2.907 


7.232 


8.837 


3.476 


-0.235 


0.802 


-0.221 






c^i 


0.001 


0.010 


0.191 


1.786 


9.852 


13.245 


9.652 


10.566 


18.5- 


-19.0 


qN 


1.172 


2.466 


5.727 


3.312 


1.524 


1.445 


0.919 


1.313 






c^i 


0.002 


0.030 


0.303 


1.617 


4.922 


5.505 


5.010 


5.045 


19.0- 


-19.5 


qN 


1.061 


1.547 


4.182 


6.443 


0.796 


-0.019 


1.628 


0.578 






cri 


0.004 


0.048 


0.312 


1.922 


6.374 


8.458 


7.330 


7.020 


19.5- 


-20.0 


qN 


0.939 


1.667 


4.892 


5.341 


2.004 


1.332 


-0.124 


0.557 






cri 


0.005 


0.125 


1.162 


4.906 


8.455 


6.342 


6.671 


6.514 


20.0- 


-20.5 


qN 


1.103 


1.815 


4.821 


1.963 


0.413 


0.279 


0.488 


0.567 






0.007 


0.174 


1.108 


3.575 


5.015 


4.066 


3.861 


3.404 




total 


Qn 


1.03 


1.55 


3.10 


2.70 


0.72 


0.25 


0.07 


0.14 






Si 


0.00 


0.05 


0.34 


1.18 


2.59 


2.12 


1.59 


1.17 






S2 


0.18 


0.39 


0.54 


1.51 


0.58 


0.26 


0.14 


0.17 



